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We characterize the macroscopic attractor of infinite populations of noisy maps subjected 
to global and strong coupling by using an expansion in order parameters. We show that for 
any noise amplitude there exists a large region of strong coupling where the macroscopic 
dynamics exhibits low-dimensional chaos embedded in a hierarchically-organized, folded, 
infinite-dimensional set. Both this structure and the dynamics occuring on it are well- 
captured by our expansion. In particular, even low-degree approximations allow to calculate 
efficiently the first macroscopic Lyapunov exponents of the full system. 

§1. Introduction 

It has now been some fifteen years since the discovery of collective behavior 
emerging out of infinite populations of incoherent nonlinear/chaotic units. In many 
instances, collective cycles have been found and studied, especially in locally-coupled, 
spatially-extended systems. Although originally discussed rather earlyP robust 
nontrivial collective chaos has been documented only in the case of globally-coupled 
populationsP^^^'®^ and the nature and dimensionality of this chaos is still a 
matter of debateP^ Indeed, many intermediate situations have been uncovered be- 
tween the infinite-dimensional chaos of weakly-coupled units and the low-dimensional 
chaos of the fully synchronized regimes occurring in the strong coupling limit. 

The effect of microscopic noise on such systems is a rather new topic of inves- 
tigationPHJ^n) Concerning collective chaos, Shibata et alP^J showed that its dimen- 
sionality can be drastically reduced by microscopic noise added to populations of 
weakly-coupled chaotic maps. Here we approach the same issue, starting from the 
fully-synchronized, strong-coupling limit. We have shown recently^' that this triv- 
ial collective chaos can be unfolded by the action of noise, allowing extra degrees of 
freedom to perturb the macroscopic dynamics. Starting from the fully-synchronized 
noiseless regime, we study in some detail the apparently low-dimensional chaotic 
regimes appearing when the noise intensity is increased. 

This paper follows a series of previous works in which, in particular, we have 
introduced a systematic order parameter expansion which approximates with in- 
creasing accuracy the macroscopic dynamics observed at strong coupling values by 
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low-dimensional effective maps for a small set of macroscopic observablesP^'^J^ 
Here we use this expansion scheme to tackle the issue of the actual dimensionality 
of the collective chaos observed. 

In the following section, we present numerical simulations illustrating the emer- 
gence of deterministic collective behavior in the infinite-size limit and briefly review 
the order parameter expansion by which we can approximate the macroscopic dy- 
namics through a set of hierarchically-organized low-dimensional systems. We then 
explore some consequences of such a representation in the case of maximal cou- 
pling. Notably, we show how different microscopic features, such as the shape of the 
noise distribution or the nonlinearity of the individual dynamics, translate into the 
collective bifurcation diagram. 

In the remaining of the paper, we address the question of the nature and extent 
of the information about the collective behavior that we can extract from our low- 
dimensional approximations. 

In Section 3, we show that the macroscopic attractor has a folded structure 
that is organized with the same hierarchy as our order parameters: the attractor 
is an infinitely-folded object whose folds on a smaller scale are captured by order 
parameters of increasing degree. 

In Section 4, we show that the hierarchy of the reduced systems also reflects 
the order in which the macroscopic degrees of freedom emerge and determine the 
collective dynamics. The reduced systems reproduce the dominant exponents of the 
Lyapunov spectrum for the population, as computed from the so-called "nonlinear 
Perron- Frobenius operator" EUEJCBJEJ i n particular, we show that, in the region 
of interest, collective chaos is characterized by only one positive Lyapunov exponent, 
while the other macroscopic degrees of freedom play the role of fast modes with neg- 
ative associated Lyapunov exponents. Nevertheless, the dimension of the attractor, 
as measured by the Lyapunov dimension, can increase when the magnitude of the 
negative eigenvalues reduces. This typically happens when the coupling is weakened. 

In the last section, we discuss the perspectives of our work and in particular the 
region of intermediate couplings and weak noise, where the collective dynamics is 
high dimensional also in the noiseless case and the order parameter expansion to low 
degree diverges. 

§2. Macroscopic dynamics 

2.1. Order parameter expansion 

We address here the population of globally-coupled noisy maps: 

Xj ^{l-K)f{x 3 )+K{f(x))+i 3 {t) j = l,...N, (2-1) 

where / : R — > R defines the individual (local) dynamics, K E [0, 1] quantifies the 
coupling of every population element to the average iterate (f{x) ) and £j(t) is an 
independent noise term, drawn according to a given distribution of finite moments 
m q . The collective dynamics is addressed, without loss of generality, in terms of the 
evolution of the simplest macroscopic variable, i.e. the average population state or 
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mean field: 

1 - 

3=1 

Numerous studies of the noiseless case have revealed that the mean field of such a 
population can display a wealth of dynamical regimes, ranging from one-dimensional 
chaos for strong coupling, when full synchronization is attained, to higher dimen- 
sional chaos in the clustering region (the dimensionality depending on the number 
of clusters) ED M S3 E3 M Ml Q ml JHJ H3 m) 

The addition of noise to this system — as well as parameter diversity — has been 
often conceived as a means of testing the genericity of the aforementioned results 
and the structural stability of the collective dynamics Another approach consists 
in fixing the parameters of the individual map and systematically vary the noise 
intensity. Changing the variance of the noise distribution has revealed that the 
stochastic terms interact with the nonlinearities of the dynamical system. For large 
populations, noise affects the high dimensional collective dynamics at low coupling 
by reducing the dimensionality of the turbulent motionP^ 

In spite of the trivial simplicity of the macroscopic dynamics in the fully syn- 
chronous regime, the modifications of the collective behavior induced by noise are 
far from straightforward. For instance, as shown by Kuramoto and Teramae, even 
weak noise can cause the population distribution moments to scale anomalously with 
respect to the noise distribution moments over a range of scalesPJE3 

In recent work, we have focused on how the addition of noise to a synchronous 
regime modifies the collective dynamicsP^ES Somewhat counterintuitively, al- 
though noise blurs the trajectory of any population element, it affects the mean field 
of infinite populations in a deterministic way. This is because, in the limit of infinite 
population size, the macroscopic dynamics is the result of the averaging of infinitely- 
many independent stochastic processes. More surprisingly, the macroscopic attractor 
still appears low-dimensional even away from the full-synchronization regime. 

As an illustrative example, we consider now a population governed by Eq. (|2-1|) . 
where 

/(x) = l- ax 2 (2-3) 

is a chaotic logistic map (a = 1.57, one-band chaos) and the noise distribution 
is taken uniform. Varying the variance of the noise distribution a 2 as a control 
parameter, a deterministic bifurcation diagram emerges when the population size N 
is increased (Fig. 

We recently introduced an order parameter expansion, which allows to approxi- 
mate, in the strong-coupling region, the deterministic collective dynamics by means 
of a low-dimensional mapP^ This macroscopic effective dynamical system is ob- 
tained by writing the evolution equation for the mean field X as coupled to other 
macroscopic variables, or order parameters: 

Q q :={e q ) g£N, (2-4) 

where e,- = Xj — X is the deviation of a population element from the mean field. In 
this case, Eq. 1)2-4(1 identifies the moments of the population distribution at a fixed 
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Fig. 1. Bifurcation diagrams of the mean field Eq. 12-21 1 for K = 1 and populations of increasing 
size TV: a) TV = 10 2 b) TV = 10 3 c) TV = 10 4 d) TV = 10 6 and reduced system. The noise intensity 
a 2 plays the role of control parameter, while the parameters of the local dynamics Eq. arc 
kept fixed. 



time instant. The infinite-dimensional map defining the evolution of all such order 
parameters is formulated in powers of (1 — K), so that a truncation to a given degree 
n provides a reduced system of dimension n which reads: 

f(X) + E n q =2 A q (X)f2 q 

(2-5) 

(i q » m q + J2 q i=1 (f){l-Kym q . i r i (X,n 2 ,...f2 n ) q = 2,...n 

where A q and Ji depend on the first IP moments of the noise distribution Of the 
remaining order parameters, those up to q = n P (P being the degree of the the map 
/) have a dynamics slaved to the the first n, while the others are constantly equal 
to the noise distribution moment of corresponding order. 

In the Appendix, we write the reduced systems up to the fourth degree when 
the logistic map Eq. 1)2-3(1 defines the individual dynamics. There, the A q and are 
explicitly derived as functions of the first order parameters. 

The population-level parameters that figure in Eqs. (|2-5|) are those of the individ- 
ual map, defining its degree of nonlinearity, the moments m q of the noise distribution, 
and the coupling constant K. The noise variance a 2 , that we typically choose as 
the bifurcation parameter, hence appears naturally as one of the parameters of the 
macroscopic description. 
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2.2. Interplay of map nonlinearities and noise moments 

One first consequence of representing the macroscopic dynamics by means of 
Eqs. (|2-5j) is that we can disentangle the effects of the local nonlinearities from those 
of the noise. Let us consider the case of maximal coupling K = 1. Here, the zeroth- 
degree representation for n — > oo is exact and takes the form: 

OO -. 

X»f(X) + V ~V*f(X)m q , (2-6) 
U 1 ] 

where T> q f is the q-th. derivative of the individual map. 

When noise is added, the uncoupled map is perturbed by a term where the 
moments of the noise distribution are weighted by the nonlinearities of the map. 
Any (polynomial) map of maximal degree of nonlinearity p will only be influenced 
by the first p moments of the noise distribution. As a consequence, in the case of 
quadratic maps such as the logistic one, the macroscopic dynamics will be affected 
in the same way by any noise distribution of given variance and is described by: 

X i-> 1 - a a 2 - a X 2 . (2-7) 

This is, of course, another logistic map which can be rescaled to the uncou- 
pled one via a simple change of variables (see Fig. ^d) for its bifurcation diagram, 
overlapping, up to finite-size effects, that of the full population). 

In order to explore how the macroscopic dynamics is affected by noises of the 
same intensity (defined as the variance a 2 ), but whose distributions have otherwise 
different features, the individual map needs to have nonlinearities of order higher 
than two. Consider, for example, the quartic map defined by a perturbation of the 
previously considered logistic equation: 

x^l-ax 2 + bx 4 , (2-8) 

where a = 1.57 and b = 0.1. For these parameter values, the quartic map is chaotic. 
The reduced system Eq. (j2-6|) reads in this case: 

X^ 1- aa 2 + bm 4 - (a-Qba 2 ) X 2 + b A 4 . (2-9) 

which is another quartic map, whose bifurcation diagram is shown in Figure [2f a) 
for uniform and Gaussian noise distributions. The order parameter expansion repro- 
duces the interaction of the single-element nonlinearities with the noise distribution 
features. 

The effect on the collective dynamics of an increase in the noise intensity allows 
us to distinguish, on a purely macroscopic basis, between microscopic stochastic pro- 
cesses with differences in the moments up to the fourth one. The distance between 
the dynamics induced by the two noise distributions are larger for larger noise in- 
tensity. Strikingly, the collective dynamics is either chaotic or stationary for strong 
noise. Figure |2] also shows that the macroscopic dynamics is bistable for interme- 
diate noise values and that the hysteresis region is slightly changed by the kind of 
noise that is applied. Similar hysteretical phenomena have been also observed in 
populations of globally coupled ODEs with independent noise,-- 1 and find here a 
simple origin within our approach. 
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For any polynomial of P- 
th degree, the macroscopic 
dynamics at maximal cou- 
pling will be exactly de- 
scribed by a polynomial of 
the same order, whose pa- 
rameters are the first 2P mo- 
ments of the noise distri- 
bution. If the uncoupled 
map is not polynomial, in 
principle all the noise dis- 
tribution moments affect the 
macroscopic dynamics and 
Eq. Q2-6[) might be a diverg- 
ing series. In some cases, 
however, the fact that m q 
scales as <r q allows to approx- 
imate the series by its trun- 
cation to a finite degree in a. 
As long as noise is sufficiently 
weak, a formal truncation to a sufficiently high degree is able to describes in detail 
the macroscopic bifurcation diagram of the full system. For instance, this applies to 
population of excitable maps, where microscopic noise gives rise to complex collective 
behavior E3J32j 

§3. Fine structure of the macroscopic attractor 

When K < 1, the zeroth-degree reduced system is no longer an exact solution 
for the collective dynamics, but Eq. (|2-5|) provides a hierarchy of maps ordered in 
powers of (1 — K). For not too small coupling values, the macroscopic dynamics still 
looks rather low-dimensional. As we have shown already in Refs. IT2^1 and !15j) . this 
hierarchy, when truncated, can account quantitatively for the collective dynamics. 
The rest of this paper is dedicated to explore in some detail what are the characteris- 
tics of the macroscopic dynamics that are captured by such low-dimensional systems. 
In particular, we will address what changes in the mean field attractor are induced 
by a reduction of the coupling strength K. 

Figure |3] displays a bifurcation diagram for fixed noise intensity and varying 
coupling strength. As for the case, illustrated in Fig. [2(d), of fixed K and varying 
a 2 , the collective dynamics seems to undergo macroscopic bifurcations among low- 
dimensional regimes. However, the scaling in the width of the chaotic bands and 
the distribution of the periodic windows suggests that the diagram might not be 
straightforwardly rescaled to that of the local scalar map. The anomalies become 
more evident when the noise intensity is weakened: for instance, for low coupling 
the macroscopic dynamics enters a region of multistability, probably associated with 
the presence of clusters in the limit of zero noise. 



x 




(T 



Fig. 2. Bifurcation diagrams of Eq. 12-91 for a local dynam- 
ics defined by Eq. 12-81 with different noise distributions: 
Gaussian (grey), uniform (black). The diagrams are in- 
distinguishable, up to finite-size effects, from those of 
the corresponding population. 
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In this Section, we study the 
structure of the macroscopic at- 
tractor and compare the simula- 
tions for the full system to those 
of the reduced systems to the first 
four degrees. As a case study, we 
choose the population of logistic 
maps whose local dynamics is de- 
fined by Eq. (|2-3|) and a uniform dis- 
tribution for the noise term. The re- 
duced systems relative to this pop- 
ulation are given in the Appendix. 

Let us start by considering the 
first return map of the mean field. 
When the coupling is maximal, the 
mean field evolves according to the 
one-dimensional map Eq. I|2-7|) . Its 
first return map hence lies onto a 
one-dimensional manifold, the parabola: 




Fig. 3. Bifurcation diagram with respect of the 
coupling strength K for the mean field of a 
population of N = 2 20 globally coupled logis- 
tic maps with added microscopic uniform noise 
of variance a 2 — 0.03. 



X(t + 1) = 1 - a a 2 -aX 2 (t). 



(3-1) 



Reducing the coupling, this parabola folds (Fig. Ufa)), which indicates that the 
macroscopic attractor is no longer strictly one-dimensional. Eq. (J2-7|) is not able to 
account for the modification of the attractor, but the folding is captured with great 
accuracy by the reduced system of second degree Eq. (jA-lj) . 




Fig. 4. a) First return map for the mean field of the population (black dots) and for the reduced 
system of second degree Eq. I|AT(I (grey dots), for K — 0.4 and a 2 — 0.001. The dashed line 
is the invariant parabola Eq. (13- 11 where the single-element dynamics is embedded, the solid 
line is the invariant parabola for the zeroth-degree approximation Eq. 12-71 1. b) The third order 
parameter f2$ versus the second Q2 for the population of logistic maps (black dots) and for the 
third degree reduced system Eq. 1A-3L The continuous line is the invariant manifold which 
embeds the dynamics of Eq. 1A-H and of Eq. 12-71 (square). 
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In spite of the fact that the reduced system of second degree reproduces qual- 
itatively and quantitatively the macroscopic dynamics of the population, this ap- 
proximation is not able to track modifications of the mean field attractor related to 
the asymmetry of the population distribution around its average. Such asymmetry 
appears as the coupling is less than maximal, and is best revealed by looking at the 
third order parameter as a function of the second order parameter Q2 (Fig. HI (b)). 
In the zeroth and second-order approximations, the dynamics of Q2 is restricted to 
the line = 0, and in the first case it reduces to the trivial fixed point (a 2 , 7713) of 
Eq. (|2-7|) . The reduced system of third degree is a three-dimensional map, reported 
in Appendix, for the mean field, the second and the third order parameter. Figure 
SI (b) shows that, contrary to the lower degree truncations, it gives rise to the same 
folding observed for the population third moment. 




'2 O 



Fig. 5. a) The fourth order parameter Q4 versus the second Q2 for the population of logistic maps 
(back dots) and for the reduced system of fourth degree Eq. l|A-4fl . The continuous line is the 
invariant manifold which embeds the dynamics of the second and third degree Eqs. (lATt and 
1A-3II and of the reduced system of zeroth degree Eq. 12-71 (square) . b) First period-doubling 
bifurcation lines for the full (dots), for reduced system of second degree Eq. and fourth degree 
Eq. l|A-4fl . The zeroth degree approximation would give vertical lines corresponding to the 
bifurcation values for K = 1. 

Finally, we project the macroscopic attractor onto the plane (J?2> ^a)-, as shown 
in Fig.^a). In the reduced system at zeroth degree, the two order parameters coin- 
cide with the noise distribution moments and thus identify with the point (<7 2 ,m4). 
In the reduced systems of both second and third degree, instead, the dynamics of 
i?4 is embedded in the invariant line defined by the conservation relation Eq. (|A-2f) . 
In the reduced system of fourth degree in (1 — K), the fourth order parameter is an 
independent variable. Figure Efa) shows that its dynamics is very close to that of 
the population and displays a similar folded structure. 

As the degree of the truncation is increased, one observes a convergence in the 
dynamics of the lower-order variables. Thus, the reduced system of fourth degree 
does not visibly improve the agreement of the lowest degree order parameters with 
the macroscopic variables of the population. As expected, the fourth degree approx- 
imation is significantly better only at low coupling values, close to the clustering 
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region. This is clear at weak noise values, when this reduced system provides a bet- 
ter approximation for the macroscopic bifurcations diagram than the lower degree 
ones (Fig. @b)). 

The distance between the full system attractor and that of the reduced system of 
lowest degree can be computed as the mean square distance of X from the parabola 
defined by Eq. l)3Tjl ). It scales as (1 — K) 2 for fixed noise intensity in a large 
region of the parameter space. 12 ' This confirms that the approximation to second 
degree indeed captures the most important deviation of the population dynamics 
from the reduced system of zeroth degree. Similarly, the distance of the fourth order 
parameter from the line Eq. (|A-2f) . valid for the approximations to second and third 
degree, scales as (1 — K) A when the coupling strength is varied. 

We conclude that our order parameter expansion describes well the fine structure 
of the macroscopic attractor. In particular, the fact that we can unveil finer and finer 
folds by adding higher order parameters suggests that the hierarchy of the folds is 
well captured by the hierarchy of the reduced system. We now turn our attention to 
the dynamics occurring on this complex object, trying to quantify the dimensionality 
of the macroscopic chaos observed and to assess to what extent it is reproduced by 
our approximations. 



§4. Lyapunov spectrum of the macroscopic dynamics 



In this section, we show that the Lyapunov exponents of the reduced system 
correspond exactly to the largest macroscopic Lyapunov exponents of the full system. 

These largest macroscopic Lyapunov exponents are most easily computed from 
the dynamics of the full system expressed in terms of the evolution of the probability 
distribution function (pdf) pt{x) of all elements at time t. This pdf is governed by 
the so-called nonlinear Perron-Frobenius operator^^^'^ 



pt+i(x) = I G{F t (y)-x} Pt (y)dy (4-1) 

where G{z} is the noise distribution and Ft the iterate of the map: 

F t (y) = (l-K) f(y)+K f p t {z)dz. 

Numerically, this is implemented easily by choosing a sufficiently fine binning of the 
support of the pdf, which is then represented by a vector. The Lyapunov exponents 
are computed from this discretized system with the method of Benettin et al.p^ 
following the evolution of a step perturbation of the pdf, according to the formula: 

1 T 

A 4 = -5>g(4), (4-2) 
t=l 

where A\ is the ratio between the Euclidean norm of the i-th component of a vector 
basis at time t — 1 and its corresponding component of the basis at time t. This 



10 



S. De Monte, F. d'Ovidio, E. Mosekilde, H. Chate 



last basis is recomputed at each time step, after a sufficiently long transient, by 
orthonormalizing the image of the basis of the tangent space at time t — 1. 

We compute the Lyapunov exponents of the reduced system by evaluating explic- 
itly the Jacobian of the system and by computing the iterates of vectors belonging 
to the tangent space. For the case of the Perron-Frobenius operator, instead, we 
iterate perturbations of small size along a trajectory. The perturbation size has 
to be chosen such that the dynamics remains close to the linear regime and a suf- 
ficient numerical precision is maintained. The perturbation size appeared not to 
be a critical parameter, except for the region close to maximal coupling K = 1. 
In this region there are directions con- 
verging very fast, so that one time 
step may be sufficient to greatly flat- 
ten the small-scale structures of a dis- 
tribution, leading to a decrease of nu- 
merical precision. This effect appears 
as a small deviation between the full 
system and the reduced one for the 
higher moments when K approaches 
1. In the following, we have checked 
that our results do not depend on the 
level of discretization of the support, 
insuring the extensivity of the spec- 
trum. 

In Fig.|f)]we show the eigenvectors 
(in the space of the pdfs) relative to 
the dominant Lyapunov exponents in 
a typical case. The increasing steep- 
ness of these objects as one goes down 
the Lyapunov spectrum signals that a 
reliable determination of many expo- 
nents is a computationally demand- 
ing task. 

In Fig. Ufa), we show the four 
largest Lyapunov exponents for the 
full system and for the reduced sys- 
tem of fourth degree as a function of 
the coupling strength K and for fixed 
noise intensity. Clearly, the reduced system captures the main features of the macro- 
scopic dynamics, i.e. the existence of only one positive Lyapunov exponent (mag- 
nified in Fig. E[b)). The quantitative agreement with the population is extremely 
good for the first two exponents. The discrepancy on the third and fourth exponent, 
evident for couplings close to the maximal one, are due to a degradation of the nu- 
merical accuracy when some directions are too strongly attractive. The quantitative 
correspondence between the population and the reduced system, especially as far 
as the third and fourth exponents are concerned, is also lost for too low values of 
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Fig. 6. Eigenvectors relative to the four largest 
Lyapunov exponents for the population of 
noisy logistic maps (K = 0.5 and a 2 — 0.03). 
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K K 

Fig. 7. a) Lyapunov exponents of the population (dots) and of the reduced system of fourth degree 
(solid lines) as a function of the coupling strength, for fixed noise intensity a 2 = 0.03. b) 
Maximal Lyapunov exponent of the population (dots) and of the reduced systems of second 
(solid line), third (dashed line) and fourth degree (dotted lines) as a function of the coupling 
strength, for fixed noise intensity a 2 — 0.03. The fourth degree truncation provides a slightly 
better approximation for the lower values of the coupling in the displayed range. 



the coupling, where the introduced approximation is likely to be insufficient to take 
every detail into account. 

In the K — ► 1 limit the maximum 
Lyapunov exponent tends to the value 
which can be derived from the exact 
scalar map Eq. (|2-6f) , while all the other 
directions are all the more attracting as 
K is close to 1. The truncations of de- 
gree greater than zero provide the suc- 
cessive exponents, each approximation 
level adding one new exponent which 
is more negative than those retrieved 
with the truncations to lower degree. 
This can be seen in Fig. |H1 which dis- 
plays the Lyapunov exponents of the 
reduced system at different truncation 
levels. 

In view of the obtained quantita- 
tive and qualitative agreement with the 
Lyapunov exponents of the population, 
we conclude that, in spite of being just 

one of the possible projections of the population dynamics on a lower dimensional 
space, the order parameter expansion indeed captures the hierarchical structure of 
the macroscopic dynamics. 
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Fig. 8. Lyapunov exponents of the reduced sys- 
tem to second (two solid lines) , third (three 
dashed lines) and fourth degree (four dotted 
lines) as a function of the coupling strength, 
for fixed noise intensity a 2 = 0.03. 
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4.1. Dimension of the macroscopic attractor 

The computation of the Lyapunov exponents also allows us to address the deli- 
cate problem of how the macroscopic attractor dimension depends on the population 
parameters K and a. The question of whether the collective behavior is truly low- 
dimensional, or instead always infinite-dimensional, has been mainly addressed in 
the context of weakly-coupled noiseless maps Shibata, Chawanya and 

Kaneko pointed out that, for weak coupling, it is sufficient to add a small amount of 
noise to reduce drastically the dimensionality of the collective motion, measured as 
the number of positive Lyapunov exponents of the spectrum of the Perron-Frobenius 
operator 

Here we approach this question starting from maximal coupling (full synchro- 
nization) rather than from zero coupling and compare the attractors of the full and 
reduced systems on the basis of their Lyapunov or Kaplan- Yorke dimension: 

D L =j+ ,, fc=1 | (4-3) 
l A i+i| 

where j is the smallest integer such that, if the Lyapunov exponents are listed 
in descending order, Yjk=x > 0. The Kaplan- Yorke conjecture, which states that 
the fractal dimension of a strange attractor is equal to appears to hold well 

for sufficiently nonsingular maps, and has been proved to be an upper estimate of 
the fractal dimension for a class of mapsP^ We thus expect to be able to use the 
Lyapunov dimension as a measure of the dimension of the macroscopic attractor. 

Although the addition of more 
and more macroscopic degrees of 
freedom improves the description of 
the mean field attractor by revealing 
smaller scale structures, as we dis- 
cussed in Section |21 this does not nec- 
essarily lead to an increase in the di- 
mension of the macroscopic attrac- 
tor. Indeed, in the region of inter- 
est here, the macroscopic attractor is 
always characterized by one positive 
Lyapunov exponent and its increase 
in dimension is associated to changes 
of the attractiveness of the transver- 
sal modes. If the coupling is suffi- 
ciently strong, the dimension of the 
attractor is less than two and is de- 
termined only by the first one of these 
modes, that associated with the sec- 
ond order parameter. 

In Fig. we show that the Lya- 
punov dimension of the reduced sys- 
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Fig. 9. Lyapunov dimension of the macroscopic 
attractor for the population (dots) and the re- 
duced system to the fourth degree (solid line) . 
The scan corresponds to the bifurcation dia- 
gram Fig. 01 for a 2 = 0.03. The values of the 
Lyapunov dimension of the reduced systems of 
lower degree nearly overlap with those of the 
fourth degree, and therefore are not displayed. 
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tern has the same dependence on the population parameters as that of the popula- 
tion. In particular, this dimension tends to one for K — ► 1, that is where we have 
analytically concluded that the macroscopic attractor is strictly one-dimensional. 
Weakening the coupling, instead, one observes an increase in the attractor dimen- 
sion, until the system exits the chaotic region through a macroscopic bifurcation, and 
the Lyapunov dimension drops to zero. Similarly, the intervals within the chaotic 
region in which the dimension of the macroscopic attractor is zero correspond to the 
windows of periodic behavior that one can identify in the bifurcation diagram shown 
in Fig. |21 

§5. Discussion 

We have studied the macroscopic attractor of large populations of noisy maps. 
We have shown that, for sufficiently strong coupling, there exists a truncation of 
the order parameter expansion which is able to reproduce the collective dynamics 
to any given scale of description. For a fixed coupling, the macroscopic dynamics is 
effectively described by a sufficiently large number n of order parameters. In this 
approximation, some of the following order parameters are slaved to the first ones, 
while all the others can be approximated by the noise distribution moments within 
the error of a macroscopic description. The hierarchy of reduced systems obtained in 
such a manner reflects the fractal structure of the mean field attractor, that can be 
described on a finer and finer scale as the degree of approximation is increased. In 
spite of its complex structure embedded into an infinite-dimensional phase, however, 
the macroscopic attractor has a finite Lyapunov dimension, which is the same as 
that of the reduced systems. 

Some further questions that can be addressed by means of the reduced systems 
are whether the macroscopic attractor attains its maximum Lyapunov dimension for 
fixed noise intensity and how this maximum value changes with a 2 . In particular, 
it would be interesting to determine whether this maximal dimension converges to a 
finite value in the limit of weak noise for any coupling strength. 

From the order parameter expansion, it is clear that the moments of the noise dis- 
tribution constitute the most relevant terms in the evolution of the order parameter 
of corresponding degree. In the limit of maximal coupling, thus, the order parameter 
of g-th degree scales "normally", that is proportionally to a q . For coupling less than 
maximal, these iterates are perturbed by the macroscopic dynamics, so that one can 
identify deviations from the normal scaling associated with the bifurcations of the 
mean field attractor. It is in principle conceivable that such bifurcations give rise to 
the "anomalous scaling" observed for intermediate coupling strengths, close to the 
point where the synchronous regime breaks down for a population of noiseless maps. 
Another possibility is that for intermediate couplings the effective dimensionality 
of the mean field dynamics diverges, so that no finite truncation is able to capture 
the scaling of the order parameters with the noise intensity. In order to answer 
this question, we are now exploring the region of low coupling/weak noise close to 
the turbulent regime, where the attractors for the low-dimensional reduced systems 
presented in this paper become unstable. 
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Appendix A 

Reduced systems up to fourth degree for logistic maps 

We explicitly derive the reduced systems for the case where individual maps 
have the quadratic form Eq. (|2-3j) . The explicit form of A q and are obtained 
following the procedure illustrated in the derivation of Eq. (J2-5jh 151 first, we write 
the equations for the first 2n moments and then we eliminate n 'slaved' degrees of 
freedom by noticing that certain linear combinations of the order parameters are 
invariant. 

As we have seen in the paragraph 12.21 Eq. (|2-5|) to the zeroth degree yields Eq. 
(t2~7j) . 

The second degree reduces system reads: 
X^> l-aI 2 -ofl 2 

Q 2 i-> o- 2 + (1 - K) 2 a 2 [m A - 6a 4 + {AX 2 -Q 2 + Qa 2 )Q 2 ] ' ' ' ' 

and the two conservation relations are valid: 

J? 3 =m 3 

Q A =m 4 + 6 a 2 {Q 2 ~ v 2 ) . (A-2) 
The third degree reduced system is the three-dimensional macroscopic map: 

'X^ 1-aX 2 -aQ 2 

< fl 2 i-> a 2 + (1 - K) 2 a 2 U 2 (A-3) 
f2 3 ^ -(l-A) 3 o 3 r 3 

where: 

E 2 =m 4 -6cr 4 + 6a 2 Q 2 - Ql + AX 2 Q 2 -AX f2 3 

E 3 =m 6 - 15 m 4 a 2 + (12 m 4 - 72 a 4 ) X 2 + (12 m 4 + 18 a 4 ) Q 2 + 72 a 2 X 2 Q 2 

-is a 2 n\ - 12 x 2 n\ + 2 n\ - 12 x n 2 n 3 + 8 x 3 n 3 . 

The following three order parameters are slaved to the first three macroscopic 
variables and they obey the relations: 

J?4 =m 4 + 6 o~ 2 (i?2 — o" 2 ) 
Q 5 =10o- 2 ^ 3 

Qq =itiq + 15 m 4 (i?2 — a 2 ) , 
while fi q = m q holds for q > 6. 
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The truncation of the order parameter expansion to the fourth degree is the 
four-dimensional macroscopic map: 



'X\-+ l-aX z -aQ 2 

r 2 i f-i is\2 Jl 



a 2 + (1 - Kf a 1 U 2 
Q 3 ^ (1 - Kf a 3 U 3 

n 4 ^ m 4 + 6a 2 (l- Kf a 2 E 2 + (1 - Kf a 4 E A 



(A-4) 



where: 



S 2 — J?4 

U 3 =f2 6 + 6X (f2 5 



n^ + Ax 2 n 2 + AXf2 3 



2Q 2 Q 3 



3 n 2 n 4 + 8 x 3 q 3 



12 x 2 (n\ 



Q 4 ) + 2Q, 



r 4 ={2 8 + 8X (3 {1% Sl 3 - 2 tt 2 fl 5 + Q 7 ) + 2AX 2 (Q 3 



-32X 3 {q 2 q 3 - n 5 ) + qqI n A + lex 4 n A 



3f2 A 



2 - 2 fi 2 f2 4 + fi 6 ) - 4Q 2 Q & 

2- 



The successive order parameters depend linearly on the first ones according to 
the relations: 

fl 5 =10 a 2 Q 3 

i?6 =m§ + (m,4 — 6 a 4 ) (i?2 — f 2 ) + Q- 2 [Q 4 — mf) 

Qj =35 777-4 ^3 

i?g =777g + 28 (?776 — 15 771-4 <J 2 ) (f2 2 — <J 2 ) + 707774 (^4 ~ "M) 

and the remaining order parameters are equal to the noise distribution moments. 
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